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co ! Abstract 

During inflation in the early universe, the Higgs field continuously acquires long-wave quantum 
" fluctuations. They accumulate to yield a non-vanishing value with an exponentially large correla- 

tion length. We study consequences of such Higgs condensations to show that, in inflation models 
^■i where the universe is reheated through gravitational particle production at the transition to the 

kination regime, they not only contribute to reheat the universe but also act as a curvaton. Un- 
fortunately, however, for parameters of the Standard Model Higgs field, this curvaton produces 
CNj " density fluctuations too large, so the inflation models followed by a long kination regime are ruled 
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I. INTRODUCTION 



The recent discovery of the new Higgs-like scalar particle by ATLAS and CMS at the 
Large Hadron Collider (LHC) [lj is likely a great step toward confirming the Standard 
Model (SM) of elementary particle physics. The discovery of the spinless particle has a 
profound implication for cosmology, whether it will actually prove to be the SM Higgs field 
or not. This is especially true for inflation in the early universe [2], which is an indispensable 
ingredient of modern cosmology to solve the horizon and flatness problems and to account 
for the origin of density fluctuations jij, since the existence of a scalar field is required to 
drive inflation. 

In this respect, the recent discovery at LHC may be directly connected with the Higgs 
inflation models [4j, where the SM Higgs field plays the role of the inflaton, the scalar field 
responsible for inflation. The SM Higgs field, however, cannot act as the inflaton adequately 
without extensions such as introducing a large and negative non-minimal coupling to the 
scalar curvature |4], a Galileon like coupling p], a kinetic coupling to the Einstein tensor 
0], a non-canonical kinetic term etc. In fact, all the five known variants of the Higgs 
inflation can be treated in a unified way in the context of the generalized G-inflation model 
(sj, which is the most general single-field inflation model with second-order field equations. 

Even if inflation is induced by a scalar field of another sector in the full theory, the Higgs 
field must exist in the SM in any case, and it plays non-negligible roles in the evolution of 
the early universe. Indeed, a renormalization group analysis shows that the self coupling 
parameter A(/i) in the Higgs potential may become negative at some large renormalization 
scale [A, depending on the values of the Higgs mass m h , top quark mass m t , and the strong 
coupling constant as{Mz)- If the Higgs field acquires an expectation value larger than 
this critical value, it would exhibit a runaway behavior to a large and negative potential 
energy which would not only hamper inflation, but would also forbid the standard big bang 
cosmology 

In December 2011, some hints of experimental signatures of the SM Higgs field were 
reported by ATLAS and CMS, with its preferred mass range 115GeV < < 131GeV by 
ATLAS and mh < 127GeV by CMS. Motivated by this result, a precise stability analysis 



of the Higgs vacua was done in [1QJ], in the range 124GeV < < 126GeV; the study 
found that for mh = 124GeV, the SM vacuum is at best metastable with a sufficiently long 
lifetime, and that the potential develops an instability at field value 10 9_14 GeV, depending 
on the values of m t and as- With a larger value of the Higgs mass, m h = 126GeV, they find 
a finite parameter space of m t and as within 2a where the SM vacuum is stable, with A(/z) 
being positive for fi at least up to the Planck scale. Even in the case where the vacuum is 
metastable, the scale where A(/z) becomes negative is pushed to a much higher scale than in 
the case with = 124GeV. Therefore, it was really fortunate for cosmology that the two 
collaborations pinned down the mass at a higher value m/, ~ 126GeV, allowing more room 
for the Higgs field value without any runaway behavior. 

In this paper, we study cosmological consequences of such a Higgs field, assuming that the 
self coupling remains positive up to a field value well above the scale of inflation. During infla- 
tion, the Higgs field acquires a nonvanishing value by accumulating quantum fluctuations, 
a process which we will call the Higgs condensation. The fluctuations are exponentially 
stretched by the subsequent inflation, and its long-wave modes become indistinguishable 
from a homogeneous mode. As a result, the universe will be filled with many exponentially 
large coherent domains, each with different field values. 
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The typical value and correlation length of the condensate can be calculated using the 



stochastic inflation method [11] . We find that for a natural value of the self coupling at the 
scale of inflation, the Higgs field acquires a fairly large value with an exponentially large 
coherent length compared to the Hubble radius during and just after inflation. Thus, when 
we calculate its effect on the reheating process, we can treat it as a homogeneous field, which 
starts coherent oscillation as its effective mass surpasses the Hubble parameter. On the other 
hand, the Higgs field is by no means fully homogeneous, so in the case where its condensate 
or its decay product significantly contributes to the total energy density of the universe, its 
fluctuation can also contribute to curvature perturbations through the curvaton mechanism 
12 1 . We show that this is indeed the case for k- inflation [13] and kinetically driven G- 



inflation [14| . in which the universe is reheated through gravitational particle production. 
Our result would also apply to the quintessential inflation [151], although it is difficult to 
realize this with a simple polynomial potential. 

The rest of the paper is organized as follows. In Sec. II, we calculate the basic properties 
of the Higgs condensation using the stochastic inflation method. In Sec. Ill, gravitational 
reheating after k-inflation is calculated using standard techniques. In Sec. IV, the contri- 
bution of the Higgs field fluctuation on curvature perturbations is discussed. 1 



II. HIGGS CONDENSATION IN INFLATIONARY COSMOLOGY 

We focus on the real neutral component of the SM Higgs field and denote it by (p(x,t). 
At the energy scale of inflation and reheating, which we assume is much larger than the 
electroweak scale, the potential of the Higgs field is well approximated by 

V(<P) = (1) 
The self copuling A = A(/z) has a logarithmic dependence on the energy scale (or field 



value). For the typical scale of inflation, we expect A(/x) = C(10 -2 ) [lOj, ll8|]. Here we assume 
that A is constant, and take its reference value as A = 0.01 hereafter. We also assume that 
the Higgs field is minimally coupled to gravity [l 



The behavior of such a scalar field during inflation has been adequately studied in [20 
using the stochastic inflation method, according to which the one-point probability distribu- 
tion function (PDF) of <p, denoted by pi(<p, t), is determined by the Fokker-Planck equation, 

^ = ^iv^ )]M } + ^ d ^$^ S K M (2 ) 



1 The possibility of the Higgs field acting as a curvaton in an inflation model based on asymptotically safe 



gravity has been discussed in 16j. A different possibility in which the Higgs field affects the density 



perturbation has recently been analyzed in fl7| , where the Higgs field is always a minor component of the 
cosmic energy density. 
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Its generic solution can be expanded as 

4n 2 V(<p)\^ * , , 



Pi(<p,t) = exp 



-An(t-to) 



3# 4 y 

Picq(v 3 ) + exp 



n=0 

4ir 2 



■A„(t-t ) 



3# 4 7 ^ 

' n=l 



(3) 
(4) 



where $n(y?) is the complete ortho normal set of eigenfunctions of the Schrodinger-type 
equation, 



1 d 2 



2 dip' 



+ W{<p) 



, , . 4vr 2 A n 



with 



and 



w%) = - K(^) 2 - , 

4vr 2 



(5) 
(6) 
(7) 



For V{ip) = ^(p we find the equilibrium one-point PDF 



32tt 2 A\ 3 



A 



exp 



2% 2 \(p 4 



3H 4 



as well as the first few eigenvalues of ((6]) as 

A = 0, Ai = 1.36859^/ ttt^H, A 2 = 4.4537 



(8) 



24tt 2 ~ 



24vr 2 



(9) 



numerically [20], which are useful in evaluating correlation functions. 

As a result, we find the following equilibrium expectation values for each quantity, 



<^ 2 > 



H 2 



Jt£M ^ ^ 0.132^= = 1.32A-^ 2 , 
2n 2 T(\)^\ v/A 



m 



cff 



A(vj 2 ) ~ 1.32 x 10" 2 A5# 2 



(^)) = ^(^) = ^^9.50xl0- 3 ^ 4 , 



(10) 

(11) 
(12) 



where we have defined A = 10 2 A ~ 1. 

Using the eigenvalues (Q, in particular the fact that A 2 is significantly larger than Ai, 
the temporal auto correlation function is well approximated by 



G(h-t 2 ) = (<p(x,ti)<p(x,t 2 )) ~ (<p 2 )e 



2\ — Ai|ti-* 2 | 



(13) 



for \ti — t 2 \ ^ H l . The correlation time, t c , which is defined by G(t c ) = |G(0), is given by 
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The spatial correlation function can also be evaluated making use of its de Sitter invari- 
ance [20]. Replacing t\ — t 2 by ^hi(Har) in flT3|) . we find 



2 A, 



G(r) = (^(x 1 ,t)^(x 2 ,t)} ~ ( V 2 )(Ha(t)r)'—, (14) 
where r = |xi — x 2 |. The spatial correlation length is defined by G(r c ) = |G(0), and reads 

a(t)r c = ^e^T ~ E_ e ™-^-\ (15) 

which is exponentially larger than the Hubble radius during inflation, as well as that at the 
end of reheating, if not larger than the comoving scale corresponding to the current horizon. 
Hence one can regard the field configuration to be practically homogeneous when we discuss 
the reheating processes. 

This means that all the massive fermions and gauge bosons in the standard model acquire 
a mass squared proportional to (y? 2 ). For quarks and leptons the resultant masses are 
substantially smaller than the Hubble parameter (except for the top quark) because they 
are suppressed by the Yukawa coupling. 

On the other hand, gauge bosons could have a mass close to the Hubble parameter. It 
breaks the conformal invariance of the vector field, and each species of gauge bosons may be 
created gravitationally with an amount close to the Higgs boson. This issue will be discussed 
elsewhere. Here we incorporate it in the effective number of modes of created particles, N, 
to be defined below. 



III. GRAVITATIONAL REHEATING AFTER K-INFLATION 

In this paper we consider consequences of the Higgs condensation discussed in the pre- 
vious section, in inflation models where the universe is reheated by gravitational particle 



production without being accompanied by the inflaton's field oscillation [13H15 



Specifically, let us consider k-inflation [13[ with its Lagrangian of the form 

C = K x {(t>)X + K 2 (<f))X 2 , X = —gTd^h (16) 

where is the inflaton and X is its canonical kinetic function. The fact that such a theory 
can realize inflation even without any potential term can most easily be seen by considering 
the situation where K\ and K 2 are constants with opposite sign. Then in the homogeneous 
and isotropic background, X has an attractor solution X = — 7^ > 0. The cosmic energy 
density, p, and the pressure, P, read 

DC 

p = IX— -£ = -P = constant, (17) 
oX 

which induces an exponential cosmic expansion with the Hubble parameter squared Hf n{ = 

2nG K\ 

3 Ki- 
ln this model, inflation terminates when Ki and K 2 both become positive. Then, the 

kinetic energy starts to redshift quickly and only the first term in the Lagrangian becomes 

relevant. That is, the universe is dominated by the kinetic energy of a free scalar field with 
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its equation of state w = P/p — 1. Then, the scale factor a(t) scales as a(t) oc i 1 / 3 , or in 
terms of conformal time, p, it asymptotically behaves as a(t) oc p 1 / 2 . The energy density 
decreases as a -6 , and the scalar curvature is given by R = —6H 2 (t). 

If we assume that the universe underwent an immediate transition from the de Sitter 
phase to the kinetic energy dominated phase (kination phase), the scale factor right after 
inflation can be written using conformal time p as 

a2{v) = Wl^¥ iv ~ Vo) + Wl^- (18) 

H inf\V0\ H mf\V0\ 

H in f here is the Hubble constant during inflation and < is the value of the conformal 
time at the end of inflation. We will take a(p ) = H ^ to be equal to 1 by adjusting r] , 
which gives 

a 2 (p) = j^(p-p ) + l- (19) 

Reheating in k-inflation is realized by gravitational particle production, due to the change 
in the definition of the vacuum. Following the standard calculation of particle creation 



2ll . l22j . we consider the creation of massless minimally coupled scalar particles. We can 



calculate the Bogolubov coefficient (3^, 

• poo 

P w = — J e- 2 ^V( V )d V , (20) 



where V(p) = ^a 2 (p)R(p), and we attain expressions for the number density and energy 



1„2 

v yii) = 

density 

n 



1 



27r 2 a 3 J 



\Pu\ 2 u 2 du, (21) 



Pr=2^l IL\ 2 ^'. (22) 







The r in n r , p r stands for the relativistic particles created at reheating. 

During inflation, (3 W is proportional to uo~ 2 at high energies, so the integral for p r diverges 
logarithmically. This divergence is due to the discontinuity of R, which can be seen by 
rewriting p r as 

i rvo rvo 

Pr = -^T- A / drn / d V2 H\ m - V2\p)V'(vi)V'(v2), (23) 

where the primes denote the derivative with respect to rj. p is some arbitrary mass inserted 
for dimensional reasons (p r is independent of the value of p). 

The above-mentioned divergence is the result of the sudden-transition approximation 
from the de Sitter phase a{rf) = —-pp-^ to f|T8|) . In order to obtain a finite value for the 
energy density, we consider the case in which the universe makes a transition from de Sitter 
space to a kinetic energy dominated universe within a timescale At] = H^xq and keeping 
the Ricci scalar continuous, following Ford (22[ who actually considered a smooth transition 
from de Sitter to a power-law expansion with a(t) oc t x l 2 instead of £ 1//3 . Since a(^o) = 1 ; we 
have Ap ~ At. We will show that for an arbitrarily small but finite At, the energy density 
after inflation can be analytically calculated. 
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Let us define 



and rewrite the energy density as 



f(H inf r)) = a 2 (77), 



H inf 



1287r 2 a 



2^4"' 



where 



dx 1 / Gte 2 ln(|xi - x 2 \)V'(xx)V'(x 2 ), 



f'f-W) 2 



V(x) 



P 



(24) 



(25) 



(26) 



(27) 



x = H m fT) is the time when V becomes sufficiently smaller than 1, at which the notion of 
particle is well-defined. 

We make the following ansatz, 



/(*) 



x 



(x < -1) 

a + a\x + a 2 x 2 + a 3 x 3 (—1 < x < x — 1) 
b (x + 61) (x — 1 < x) 



(28) 



and require /(x), f'(x), f"(x) to be continuous at a; = — 1 and 2 = Xo — 1 (which makes the 
Ricci scalar continuous throughout the period under consideration). If the transition from 
Ki((j))K2((j)) — constant < to Ki((f))K 2 (<p) > occurs in a time scale much less than the 
Hubble time, numerical calculations show that the transition to the kination regime with 
w — 1 also takes place well within the Hubble time 14]. Thus we can take xq < 1, or 
At<Hrl 



The coefficients are determined as (up to 0((xq)°)), 



1 

a = 6 , 

x 

1 

03 = , 

x 



which leads to 



a% = 8 
bo = 2, 

V'=- 



x 



6 

Xq 



a 2 = 3 - 

3 

2' 



x 



61 



(29) 
(30) 

(31) 



We can see that 60 an d 61 coincide with the expression for the scale factor f fTB"]) . which 
justifies our ansatz. 

We can now calculate I, which becomes 

J~-361nx , (32) 

so the energy density at the beginning of the kinetic energy dominated era is 

1 



Pr 



UU ™ f A In 
32vr 2 a 4 



H mi At 



(33) 
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This value is % times larger than the value for the radiation dominated universe, derived in 



221 ] . We will ignore the dependence on At which is logarithmic, and take \n(l/(HAt)) ~ 1. 
If there are N modes of species created this way, (1331) should be multiplied by N. 

Now we go on to calculate the reheating temperature. If we neglect the effect of Higgs 
condensation, the reheating temperature is determined by the condition p^ = p r , where p^ 
is the energy density of the inflaton field, 

P4> = 3M 2 G H^a- e . (34) 

Here Mq is the reduced Planck mass. Then, the scale factor clr at the time of reheating 
satisfies 

_ 32^ 

R " ' (35) 

so the energy density of the relativistic particles is 

pAr ~ (32n^ G (36) 

In terms of the reheating temperature, the radiation energy density is expressed as ~T Rl 

30 

where g* is the effective number of relativistic species of particles, so 

T R = # - 3.9 x 10^! (-*.) k (JL-) GeV. (37) 

(327r 2 )f WdJ m g V106.75/ V0.01/ V ; 

Here r denotes the tensor-to-scalar ratio which is given by 

r = 0.0lf H ™ V. (38) 

V2.4 x 10 13 GeVy V ; 

The above analysis applies in the case in which the Higgs condensation does not contribute 
to the cosmic energy density. Now we consider how the Higgs condensation affects the 
reheating temperature in this model. After inflation, the Higgs field remains constant until 
the Hubble parameter decreases below m e s ~ 0.115A 1//4 i^i n f (see eq. (TTTi) ). When H becomes 
smaller than m e g, the Higgs field starts oscillating. At the beginning of the oscillation phase, 
the energy density of the inflaton, relativistic matter, and Higgs field can be expressed as 

p^ = 3M G m 2 e{i ~ 3.96Aiil4if? f , (39) 

Pcond = ^ 9.50 x 10- 3 H? n( , (41) 

respectively. Thus p C ond can be larger than p r at this moment. Since both p r and p con d scale 



as a~ after the oscillation begins 23] , the universe can be predominantly reheated through 
the decay product of the Higgs condensation. In such a case, the reheating temperature can 
be estimated from the equality p^ = p con d, because the Higgs condensation dissipates its 
energy to radiation rapidly, once the field oscillation commences. As a result, we find 

T R — 1.8 x 10 7 ( — ^— V 1 ( —) GeV. (42) 
R V106.75/ V0.01/ V ; 
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IV. FLUCTUATION OF THE HIGGS FIELD AFTER INFLATION 



So far we have treated the Higgs field as a practically homogeneous condensate after 
inflation, because it has an exponentially large correlation length compared to the horizon 
scale at the end of inflation. In reality, however, the Higgs field acquires long-wave quantum 
fluctuations whose observational consequence should be clarified. 

To do this, we first calculate its long-wave power spectrum from the spatial correlation 
function (TI4T) . A power-law correlation function can be obtained from a power-law power 
spectrum P{k) = \(pk\ 2 = Ak n as, 

G(r) = J P(k)e^^ = Ar* 3 ^±^ sin [(„ + 2)|] , (43) 

2Ai 

which, for the case we are considering, means n = —3 H and 

Hint 



An 2 A 1 / 3 r(|) _^i„2- 



2A, 



^^V^Afi ^^"- < 44) 



respectively. In deriving the latter, we have adopted an approximation that -n- 1 - is much 

^inf 

smaller than unity. 
As a result we find 



PM = W , = i*4(5) y «|, («) 

that is, its power-spectrum is very close to that of a massless minimally coupled field. 

We can also calculate the power spectrum of the energy density fluctuation of the Higgs 
condensate after inflation. Since it is frozen while m e fr < H m f, we can estimate it only from 
the potential energy fluctuation. Since the potential energy is also governed by the same 
PDF as used above, its large-scale correlation function is given by 

= a » S / «»Mfo(0) \ a 1 [{v(r)vm _ {vm 

\ Ph Ph I \V) 

' [(V 2 ) - (V) 2 ] e~ Al ** = A(Har)- 2 -^ 



(V) 

Mk)e-^ (46) 



where = 4ln(ifar). 

The power spectrum Pv(k) then reads Py(k) ~ 14a/A/c -3 , so the amplitude of density 
fluctuations on scale r = ^ is given by 

V v (k) = -—,P v (k) = 0.71 VA = 0.071A2 (47) 

[ZTT) 6 

Thus the Higgs field acquires almost scale-invariant fractional density fluctuations with 
amplitude ~ 0.1. 
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Since this or its decay products will make a significant contribution to the energy density 
of the universe, if not dominant, after the inflaton's kinetic energy density has dissipated 
away according to oc a -6 , the Higgs field unexpectedly acts as a curvaton EH. Its 
contribution to the final curvature perturbation is schematically given as [IT 



K 



Ph 

Ptot 



8p h ,SH 



Ph 



Ph 

Ptot 



(48) 



where H osc is the Hubble parameter at the onset of Higgs oscillation (= m eff ), 5H osc is its 
fluctuation, and c and d are quantities of order unity. As a result, it induces curvature 
perturbations too large with the phenomenologically inferred value of A, since at the 
onset of oscillation is 0(1). 

Thus the decay product of the Higgs condensation cannot constitute the dominant com- 
ponent of the universe. In other words, some part of the inflaton's energy must be directly 
transferred to radiation without too long a period of the kination regime. 



V. DISCUSSION 

We have considered cosmology of the SM Higgs field based on the recent discovery of a 
Higgs like scalar particle with a mass ~ 126GeV. Assuming that its self coupling X(p) 
is positive and O(10~ 2 ) at the scale of inflation, we have analyzed its spatial configuration 
using the stochastic inflation method. We have found that for an experimentally inferred 
value of the self coupling A(p) ~ 0.01 at the scale of inflation, the Higgs condensation 
due to long-wave quantum fluctuations acquired during inflation suffers from large density 
perturbations. 

In the conventional potential-driven inflation, slow-roll inflation is followed by field oscil- 
lation, whose energy dissipates in proportion to a~ 3 until reheating takes place. Hence the 
energy density of the Higgs condensation, which is much smaller than the inflaton's energy 
density with p con d/p^ — H^/Mq < 10~ 5 , never contributes to the total energy density 
appreciably, since its coherent oscillation dissipates in proportion to a -4 . 

In k-inflation or kinetically driven G-inflation as well as in quintessential inflation, re- 
heating occurs through gravitational particle production and the inflaton energy dissipates 
as a" 6 . As a result, even if the oscillation of (p is governed by the quartic potential and 
its energy density scales as a -4 , it will eventually contribute to the total energy density 
significantly and act as a curvaton. This is an unusual case, since a scalar field with such a 
quartic potential usually does not play the role of a curvaton as explained above (see also 



[24j . |25| for other effects of self interaction). 

Thus, k-, G-, and quintessential inflation with long kination regimes are ruled out due 
to the SM Higgs field. It has been argued that the imprint of the kination regime can 
be observed as an enhancement of the high-frequency part of the stochastic gravitational 
wave background energy density produced by inflation (26|. Our analysis shows that such 
a possibility is inconsistent with the observed amplitude of curvature perturbations, due to 
the extra contribution of the Higgs field. 

There is, however, a simple remedy. If radiation particles are directly created from the 
inflaton due to some direct coupling, then reheating would be more efficient. Another 
possibility is to extend the structure of the Higgs field. For example, if one introduces 
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supersymmetry, the condensate of the H u Hd flat direction may dominate the energy density 
of the post-inflationary universe eventually. Then it not only works as a curvaton but also 
creates practically all the radiation components 27|. In either case, a long kination regime 



would not be realized after inflation and the enhancement of the high frequency gravitational 
wave background would be impossible. Thus the observational significance of our analysis 
mentioned above does not change even with these solutions. 



Acknowledgments 

We would like to thank Kohei Kamada for useful communications. This work was par- 
tially supported by JSPS Grant-in- Aid for Scientific Research No. 23340058 (J.Y.), and the 
Grant-in-Aid for Scientific Research on Innovative Areas No. 21111006 (J.Y.). 



[1] ATLAS Collaboration |arXiv:12Q7.72T4l [hep-ex]]; CMS Collaboration [arXiyil207\7235 [hep- 
ex]] 

[2] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980); K. Sato, Mon. Not. Roy. Astron. Soc. 195, 
467 (1981); A. H. Girth, Phys. Rev. D 23, 347 (1981) 

[3] V. F. Mukhanov, G. V. Chibisov, JETP Lett. 33, 532-535 (1981); S. W. Hawking, Phys. Lett. 
B 115, 295 (1982); A. A. Starobinsky, Phys. Lett. B 117, 175 (1982); A. H. Guth and S-Y. Pi, 
Phys. Rev. Lett. 49, 1110 (1982) 

[4] J. L. Cervantes-Cota and H. Dehnen, Nucl. Phys. B 442, 391 (1995) [arXiv:astro-pli/9505069l ; 
F. L. Bezrukov, M. Shaposhnikov, Phys. Lett. B 659, 703 (2008) [ arXiv:0710.3755[ hep- 
th]]; A. O. Barvinsky, A. Y. Kamenshchik and A. A. Starobinsky, JCAP 0811, 021 (2008) 
|arXiv:0809.2i~04l [hep-ph]]: F. Bezrukov, D. Gorbunov and M. Shaposhnikov, JCAP 0906, 
029 (2009) [a rXiv:0812.3622l [hep-ph]]; 

[5] K. Kamada, T. Kobayashi, M. Yamaguchi, J. Yokoyama, Phys. Rev. D 83, 083515 (2011) 
|arXiv: 1012.42381 [astro-ph.CO]] 

[6] C. Germani and A. Kehagias, Phys. Rev. Lett. 105, 011302 (2010) [ arXiv:1003.2 635 [hep- 
ph]]; C. Germani and A. Kehagias, JCAP 1005, 019 (2010) [Erratum-ibid. 1006, E01 (2010)] 
[arXiv: 1003.42851 [astro-ph.CO]]; C. Germani and Y. Watanabe, JCAP 1107, 031 (2011) 
[Addendum-ibid. 1107, A01 (2011)] |arXiv: 1106. 050211 [astro-ph.CO]] 

[7] K. Nakayama, F. Takahashi, JCAP 1102, 010 (2011) [ arXiv: 1008.44571 [hep-ph]] 



11 



[8] T. Kobayashi, M. Yamaguchi, J. Yokoyama, Prog. Theor. Phys. 126, 511 (2011) 

|arXiv:1105.5723 [hep-th]] 
[9] J. R. Espinosa, G. F. Giudice and A. Riotto, JCAP 0805, 002 (2008) [ arXiv:0710.2484l [hep- 

ph]] 

[10] J. Elias-Miro, J. R. Espinosa, G. F. Giudice, G. Isidori, A. Riotto and A. Strumia, Phys. Lett. 
B 709, 222 (2012) |arXiv:1112.3022l heD-Dh]] 

[11] A. A. Starobinsky, in Fundamental Interactions, edited by V. N. Ponomarev (MGPI Press, 
Moscow 1984), p. 55; A. A. Starobinsky, in Field Theory, Quantum Gravity and Strings, edited 
by H. J. de Vega and N. Sanchez, Lecture Notes in Physics 246, 107 (1986) 

[12] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) |arXiv:hep-ph/0109214| ; D. H. Lyth 
and D. Wands, Phys. Lett. B 524, 5 (2002) [arXiv:hep-ph /0110002] ; T. Moroi and T. Taka- 
hashi, Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)] [arXiv:hep-ph/0110 096| 

[13] C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. B 458, 209 (1999) 
|arXiv:hep-th/9904075l 

[14] T. Kobayashi, M. Yamaguchi, J. Yokoyama, Phys. Rev. Lett. 105, 231302 (2010) 

|arXiv:1008.0603 [hep-th]] 
[15] P. J. E. Peebles, A. Vilenkin, Phys. Rev. D 59 063505 (1999) |arXiv: astro-ph / 98 1 0509] 
[16] Y. -F. Cai and D. A. Easson, larXiv:1202.1285l [hep-th] 
[17] A. De Simone, A. Riotto [arXiv: 1208. 13441 [hep-ph]] 

[18] J.A. Casas, J.R. Espinosa, M. Quiros, Phys. Lett. B 342, 171 (1995); J.A. Casas, J.R. Es- 
pinosa, M. Quiros, Phys. Lett. B 382, 374 (1996); G. Isidori, G. Ridolfi, A. Strumia, Nucl. 
Phys. B 609, 387 (2001) [arXiv:hep-ph/0104016] ; G. Isidori, V.S. Rychkov, A. Strumia, N. 
Tetradis, Phys. Rev. D 77, 025034 (2008) |arXiv:0712.0242l [hep-ph]] 

[19] B. L. Spokoiny, Phys. Lett. B 147, 39 (1984); T. Futamase and K. -i. Maeda, Phys. Rev. D 
39, 399 (1989); D. S. Salopek, J. R. Bond and J. M. Bardeen, Phys. Rev. D 40, 1753 (1989); 
R. Fakir and W. G. Unruh, Phys. Rev. D 41, 1783 (1990) 

[20] A. A. Starobinsky and J. Yokoyama, Phys. Rev. D 50, 6357 (1994) | arXiv:astro-ph/9407016~] 

[21] L. Parker, Phys. Rev. 183, 1057 (1969); Y. B. Zel'dovich, A. A. Starobinsky, Zh. Eksp. Toer. 
Fiz. 61, 2161 (1971) [Sov. Phys. JETP 34, 1159 (1972)]; Y. B. Zel'dovich, A. A. Starobinsky, 
Pis'ma Zh. Eksp. Toer. Fiz. 26, 373 (1977) [JETP Lett. 26, 252 (1977)] 

[22] L. H. Ford, Phys. Rev. D 35, 2955 (1987) 



12 



[23] M. S. Turner, Phys. Rev. D 28, 1243 (1983) 

[24] M. Kawasaki, T. Kobayashi, F. Takahashi, Phys. Rev. D 84 123506 (2011) [ arXiv:1107.60TT1 
[astro-ph.CO]] 

[25] K. Enqvist, S. Nurmi, G. Rigopoulos, O. Taanila and T. Takahashi, JCAP 0911, 003 
(2009) [arXiv:0906.3l26l [astro-ph.CO]]: K. Enqvist, Prog. Theor. Phys. Suppl. 190, 62 (2011) 
|arXiv:1012.17lT1 [astro-ph.CO]] 

[26] H. Tashiro, T. Chiba and M. Sasaki, Class. Quant. Grav. 21, 1761 (2004) 
[arXiv:gr-qc/ 0307068| 

[27] K. Enqvist, S. Kasuya, A. Mazumdar, Phys. Rev. Lett. 93, 061301 (2004) 
|arXiv:hep-ph/0311224] 



13 



